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Using the Good and Walker picture, we derive a simple formula for diffractive 
dissociation that can apply to recent data collected at HERA in the low Q 2 regime. 

Deep inelastic scattering of a photon of virtuality Q off a given target can be 
viewed as the photon splitting into a quark- antiquark pair of size r ~ 1/Q, and 
subsequently scattering off the target through a further quantum fluctuation. This 
interpretation is valid in the leading logarithmic approximation of QCD when the 
total rapidity Y is very large, and in a frame in which Y is mostly carried by the 
projectile photon. The rapidity of the target in that frame is Y <C Y. 

The initial qq pair is a color dipole. The higher Fock states built up from QCD 
radiation can also be interpreted as collections of dipoles [3]. That holds in the 
large N c limit in which a gluon is essentially a zero-size qq state, as far as color is 
concerned. 

There are typically two classes of events observed in experiments: either the 
target breaks up and the whole rapidity range is filled with decay products, or it 
remains essentially intact and an angular region empty of final state particles is 
seen in the detector. The latter events are called diffractive, and may represent up 
to 20% of all events. Among them, the ones that exhibit a final state made of the qq 
pair together with additional gluonic radiation in the photon fragmentation region 
are called dissociative. They are characterized experimentally by a large invariant 
mass Mx ^> Q of the decay products of the photon. 

The main goal of the work on which we report here pQ is to provide a simple 
derivation of the dissociative cross section valid also in the saturation regime, when 
effects due to unitarity corrections become sizeable, and to compare it to recent 
HERA data. We also argue that diffractive and total cross sections have the same 
rapidity dependence, and this comes about very naturally in our framework. 



1 High energy behavior of deep-inelastic scattering observables 

From a simple quantum mechanical calculation, one obtains that the total, elastic 
and dissociative cross sections for the scattering of a color dipole of size r at fixed 
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impact parameter b are given by 

=2 (1 - (S{r, b)) Y - Yo ) , §f = (1 - S(r, b))l_ Yo 

^= {S *(r,b))Y-Y -{S(r,Y)f Y _ Yo . 

The diffractive cross section is the sum of the elastic and dissociative ones. S(r, b) 
is the S'-matrix element for the elastic scattering of a fixed partonic configuration 
of the dipole of size r at impact parameter b. It is a random variable, whose 
distribution is related to the probability distribution of a given Fock state of the 
initial qq pair. 

The average (.)y-y is taken over all possible partonic configurations of the 
initial dipole after evolution over a rapidity range Y — Yq. The last formula above 
is due to Good and Walker [2J, and is an identity between the dissociative cross 
section and the variance of the S'-matrix for partonic states. The total cross section 
instead is related to the expectation value of the latter. 

The evolution law for S with Y can be readily derived. Within the rapidity 
interval dy, a dipole of size r may split into two dipoles of size z, r — z by emitting 
a gluon. That occurs with a probability proportional to dy, say Xdy. Let us denote 
by p(z, r) the distribution of the sizes z, r — z of the emitted dipoles. Then 

!S(r) with probability 1 — Xdy 
S(z)S(r - z) with probability Xdy , 

z being distributed according to p(z, r)d 2 z. That evolution law implies a recursion 
relation for (S(r)) Y : 

(S(r)) Y+dy = (1 - Xdy)(S(r)) Y + Xdy J d 2 z p(z , r) (S (z) S (r - z)) Y . (3) 

The limit dy — > yields 



d Y (S(r)) Y = —J d 2 z z2{r _ z)2 ((S(r)) Y - (S(z)S(r - z)) Y ) . (4) 
We have used the known QCD result for the dipole splitting probability 

OL T 2 

dy X x p(z, r) d 2 z = dy- tt-, r^d 2 z , (5) 

Z7T z z (r — z) z 

where ot = a s N c /n. The latter equation is not closed: it involves the correlator 
(S(z)S(r — z))y ■ It turns out to be the first equation of an infinite hierarchy named 
after Balitskh. A mean field approximation (S(z)S(r— z)) ~ (S(z))(S(r — z)) casts 
Eq.QJ into a closed form, known as the Balitskh-Kovchegov (BK) equation 0|. 

Note that the fact that S is a scattering matrix element plays no special role in 
the above derivation. The only assumption is that the created dipoles are indepen- 
dent in such a way that their interaction factorizes. Thus we see that the quantity 
(S 2 (r)) Y obeys exactly the same evolution equation @, with S replaced by S 2 . 

The BK equation falls in the universality class of the Fisher-Kolmogorov equa- 
tion, and was shown |S] to admit traveling wave solutions at large Y, that decay 
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like 1 - (S{r)) Y ~ 1 - (S 2 (r))Y ~ {rQ s {Y)) 2l ° for r « 1/Q 2 S (Y) (70 ~ 0.63 is 
fixed). Q S (Y) is the saturation scale. From Eq.Q, it results that the ratio 

Vdif, = l-2(S(r)) Y _ Yo + (S 2 (r)) Y _ Yo 

2(1 -{S(r)) Y -Y ) ' W 

which involves quantities that all obey the same equation, is independent of Y. A 
similar result was obtained by a direct calculation in Ref.[B]. 

The independence of that ratio with respect to the center-of-mass energy was 
seen experimentally in the HERA data. However, it seems to hold in rather narrow 
bins of the invariant diffracted mass. Here, we have shown that the total diffractive 
cross section (with a minimum rapidity gap Yq) is independent of Y, which is a 
weaker result: it just means that at very high energy, deep inelastic scattering off 
a proton and off a Pomeron have the same energy dependence. 

2 High mass diffraction 

We compute the diffractive cross section to order a s , allowing for at most one 
gluon in the diffractive system. Beyond that one-fluctuation, the 5-matrix is ap- 
proximated by its average S(z) ~ (S(z)). Replacing Eq.(@J for (S) and (S* 2 ) resp. 
into Eq.JTJ and setting Yq = Yp the size of the rapidity gap, one gets |7I1| 

dtJdiss 1 j ■ 

— / a z- 



d?bd\og{l/x g ) 2tt J z 2 (r-z) 2 

x {S(Y r ,z,b- {r-z)/2)S(Y P ,r - z,b-z/2) -S{Y P ,r,b)) 2 (7) 

where x g is the longitudinal momentum fraction carried by the diffracted gluon. It 
is related to the rapidity variables through log(l/sc 9 ) = Y ~ Yp, and to the diffracted 
mass by Mx = Q/^/x^. S is the dipole elastic S*-matrix element in the mean field 
approximation. 

To arrive at HERA phenomenology, one has to convolute Eq. Q with the distri- 
bution of dipoles of size r in the virtual photon, given by the squared wave function 
\ip(z q ,r)\ 2 for the splitting 7* — ► qq (x g is the longitudinal momentum fraction 
carried by the quark). We also assume the Golec-Biernat-Wiisthoff (GBW) model 
[S] for S, which in particular has no 6-dependence. After some algebra, one gets 

dM x ~ tt 2 MxJ d r J dz ^^r)\ 

r 2 

d 2 z— -(S(Y P ,z)S(Y P ,r~ z)~S(Y P ,r)) 2 . (8) 

z z (r — zy 

cto is a parameter, and two more parameters enter S. All of them are completely 
determined by the GBW fit to the F 2 structure function [H]. 

The agreement with the data so far published by the ZEUS collaboration [§] is 
good, see Fig.l. We expect that much could be learned on the saturation regime 
of QCD from observables such as da diss /d log Q 2 if that could be measured in the 
same kinematical regime ^U]. We refer the reader to Ref.pQ for further details. 
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Figure 1. Comparison of the model with the data. Dashed line: Eq.JgJ. Dotted line: elastic 
component, taken from GBW. Full line: total. Data points are from the ZEUS coll. 
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